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ABSTRACT 


This  thesis  is  concerned  with  the  degree  of  exactness 
of  Tchebycheff  and  modified  Tchebycheff  quadrature  formulae 
on  the  semi-infinite  and  infinite  intervals. 

The  first  Chapter  is  a  historical  survey,  giving  an 
outline  of  the  development  of  the  theory  of  Tchebycheff 
quadrature  and  indicating  some  of  the  areas  of  recent 
research.  In  Chapter  2  we  examine  a  modified  Tchebycheff 
quadrature  formula  with  n  nodes  on  the  interval  (-00  ,  00  ) 
and  show  that  the  degree  of  exactness  of  this  formula  is 
of  order  log  n.  The  result  is  an  extension  of  the  results 
of  Erdos  and  Sharma,  and  Meir  and  Sharma  to  the  infinite 
interval,  although  the  modified  Tchebycheff  formula  is  less 
general  than  theirs. 

Chapter  3  deals  with  Tchebycheff  quadrature  on  the 

infinite  and  semi-infinite  intervals  with  weight  function 
2 

e  x  |  x  |  a,  a  >  -  \  for  (-°°,°°)  and  the  weight  e"~xxa,  a  >  -  1 
for  (0,o°).  These  results  are  qualitative  extensions  of  the 
results  of  Salzer  and  Gatteschi.  In  Chapter  4,  a  partial 
result  is  given  concerning  the  degree  of  exactness  of  a 
modified  Tchebycheff  quadrature  formula  on  the  interval 
(0,°°)  and  a  conjecture  is  formulated. 
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CHAPTER  I 


INTRODUCTION  -  A  HISTORICAL  SURVEY 

1.1  GAUSS  QUADRATURE  FORMULA 

The  problem  of  quadrature  or  of  finding  the  area 
inside  a  given  closed  curve  has  attracted  the  attention 
of  mathematicians  since  the  time  of  Archimedes,  when  he 
found  the  area  of  a  circle  by  subdividing  it  into  small 
triangles.  The  essential  step  in  almost  all  quadrature 
formulae,  at  least  those  of  interpolatory  type,  seems  to 
be  to  replace  a  given  function  by  a  simpler  function. 

Thus  in  the  trapezoidal  formula,  we  replace  a  given 
function  by  a  broken  line  which  interpolates  the  function 
in  a  certain  number  of  points,  while  in  Simpson’s  formula, 
a  parabolic  arc  replaces  the  given  function  over  certain 
intervals.  The  same  idea  underlies  the  general  Newton- 
Cotes  formula,  where  the  function  is  replaced  by  a  suitable 
polynomial  which  interpolates  the  function  in  equidistant 
points . 

Among  the  plethora  of  quadrature  formulae  available 
today,  the  Gauss  formula  has  an  importance  of  its  own. 

This  formula  is  also  of  an  interpolatory  character,  but 
here  the  nodes  of  interpolation  are  not  fixed  in  advance 
as  is  done  in  Simpson’s  or  the  Newton-Cotes  formula.  If 
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we  construct,  for  any  continuous  function  f(x),  the 
Lagrange  interpolation  polynomial 


L(x)  =  l  f(x  )  i  (x) 
k=l  K  K 

where 

*k(x)  =  - -  ;  x(x)  =  n  (x-x  ) 

X»(xk)(x-xk)  k=l 


and  x^  <  <•••<  x^  are  arbitrarily  assigned  nodes,  we 

have  f(x)  =  L(x)  +  Rn(x)  where  Rn(x)  vanishes  whenever 
f(x)eii  the  class  of  all  polynomials  of  degree  <_  n-1. 

Indeed 


R  (x)  = 
n 


X  (x) 


n ! 


fn(5) 


X1  < 


X 


n 


Hence,  whenever  f(x)en  ,  ,  we  have 
9  n-1 9 


In  1 

(1.1.1)  /  f ( x )  dx  =  l  A,  f ( x ,  )  ,  A  =  /  i,  (x)  dx  . 

-1  k=l  K  -1 


By  appropriately  choosing  the  n  nodes  of  interpolation, 
it  is  possible  to  increase  the  class  of  polynomials  for 
which  (1.1.1)  is  exact.  The  optimum  situation  is  achieved 
when  the  x,  ’s  are  taken  to  be  the  zeros  of  P  (x),  the 
Legendre  polynomial  of  degree  n.  This  gives  the  Gauss 
quadrature  formula: 
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1  n 

(1.1.2)  /  f(x)  dx  =  l  A  f(x,  )  ,  A  >  0 

-1  k=l  K  K  K 

valid  for  all  polynomials  of  degree  <_  2n-l. 

Closely  related  to  the  Gauss  quadrature  formula  are 
those  due  to  Radau  and  Lobatto.  If  the  end  point  -1,  of 
the  interval  of  integration,  is  an  assigned  node,  we  then 
have  the  Radau  quadrature  formula.  Assigning  both  the  end 
points  -1  and  1  as  nodes,  we  obtain  the  Lobatto  quadrature 
formula. 

For  a  more  exhaustive  treatment  of  the  theory  of 
approximate  integration,  we  refer  to  Krylov  [9]  or 
Natanson  [ 1 4 ] . 


1.2  GENERALIZATIONS  OF  GAUSS  QUADRATURE 

Turan  [21]  has  generalized  the  Gauss  quadrature 
formula  by  considering  derivatives  of  the  function  in  the 
quadrature  formula.  The  existence  of  a  formula  of  the 
type : 


(1.2.1) 


1 

/  f ( x )  dx 

-1 


n 


l 

k=l 


f(xk) 


+ 


n 


l 

k=l 


f ' 


x 


(1) 

k 


follows  immediately  from  the  interpolation  formula 


f  (x)  = 


n 

l 

k=l 


f(xk) 


*k,0(x 


)  + 


n 

l 

k=l 


(xk) 


Vi(x) 
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valid  for  all  polynomials  of  degree  <_  2n-l,  where,  using 
the  notation  of  the  Lagrange  interpolation  polynomial. 


a 


k ,  0 


(x) 


x"(xk) 

FTx^ 7  (x_xk) 


(x) 


*k  =  ( x“xk ^  £k^x^  • 


For  any  choice  of  x^,  it  is  only  possible  to  have  (1.2.1) 
valid  for  polynomials  of  degree  _<  2n-l. 

Turan  has  given  a  generalization  of  formula  (1.2.1) 
in  the  form 


n 


/  f ( x )  dx  =  l  [f(xk)x£0)  +  f' (x,,)x,(1)  +•••+  (x,jx,(. 


m 


-1 


k=l 


‘k  k 


k7  k 


valid  for  all  polynomials  of  degree  ^(m^+m^t • • • +m  -1 ) . 

If  ~  *  *  ’  =  mn  =  k ,  Turan's  formula  is  valid  for 

all  polynomials  of  degree  <_  kn-1  and  may  be  made  valid 
for  all  polynomials  of  degree  <_(k+l)n-l,  when  k  is  odd, 
by  a  suitable  choice  of  the  nodes. 

Considering  this  type  of  quadrature  formula  with 
weight  function  (1-x2)-^  ,  Kis  [8]  has  given,  for  k=3, 
the  explicit  form  of  the  quadrature  formula  considered  by 
Turan  which  is  valid  for  all  polynomials  of  degree  <_  4n-l 


1.5  '  • 

.  (x)  .1  („X-X)  =  X  £<}i 

geib  lo  alalmon 

(VTOl  9rf^  fit 
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1 


/ 

-1  Jl-X2 


H 

n 


n 


r  ,  2k- 1  s 

L  g(cos  tt) 


k=l 


2n 


1  V  „  2k- 1  ,,  2k- 1 

^  k-1  ~2tr*S '(cos  ~2n~ 


7T 


,  1  v  _ j „2  2k- 1  2k- 1  N 

+  — ?  L  sln  — p-r—  it  g  ( cos  -p—  tt) 

4n2  k=l  2n  2n 


This  is  a  generalization  of  the  Hermite  formula  of  quadra¬ 
ture  ( formula  (1.4.1)). 


1.3  ERROR  ESTIMATES 

If  we  are  applying  a  quadrature  formula,  it  is  of 
interest  to  have  some  estimate  of  the  error  involved.  For 
an  arbitrary  quadrature  formula 

b  m-1 

/  f(x)  dx  =  L(f)  5  l  A.  f(x.  ) 
a  k=0  K  K 

which  is  exact  for  polynomials  of  degree  <_  r-1,  a  natural 
question  arises  regarding  the  error  when  f(x)  belongs  to 
a  class  of  functions  bigger  than  A  discussion  of  the 

various  aspects  of  the  error  estimates  for  certain  classes 
of  functions  and  the  choice  of  nodes  which  would  minimize 
the  -error  in  some  norm  would  lead  us  too  far  afield.  We 
refer  here  to  the  works  of  Nikol'skii  [15]  and  Sard  [17]. 
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1.4  TCHEBYCHEFF  QUADRATURE 

In  1867,  Hermite  gave  the  formula 


(1.4.1) 


2k- 1  s 

“TrT  *> 


which  is  exact  for  all  polynomials  of  degree  <_  2n-l  and 
so  is  a  Gauss-type  formula.  Tchebycheff  noted  the  advantage 
of  this  type  of  formula  in  numerical  computation  and  in  1873 
posed  the  following  problem: 

Find  a  real  A  and  real  nodes  x^(l<_k_<n)*  in  such  a 
way  that 


1 


n 


/  f  ( x )  w  ( x )  dx  =  A  l  f  ( x,  ) 
-1  k=l  K 


(1.4.2) 


is  exact  for  all  polynomials  f(x)  of  degree  <_  n.  It  is 
only  necessary  to  consider  the  interval  of  integration  to  be 
[-1,1]  since  any  finite  interval  [a,b]  can  be  reduced  to 


2xr(a+b)  -*■  x .  a 


[-1,1]  by  the  linear  transformation 


Tchebycheff  quadrature  formula  is  said  to  be  possible  for 
n  with  weight  function  w(x)  if  there  exist  numbers 
x^. (k=l , 2 , .  .  .  ,n)  lying  interior  to  [-1,1]  such  that  (1.4.2) 
is  true  for  every  polynomial  of  degree  <_  n. 


*  It  should  be  emphasized  that  x^  =  x^  and  A  e  A 
depend  on  n,  but  for  typographical  reasons  we  shall 
here  and  throughout  the  thesis  omit  the  superscript  n 
With  the  understanding  that  this  dependence  on  n  is 
implied  in  the  formulation  of  the  problems  and  the  state¬ 
ments  of  the  results. 
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In  connection  with  the  formula  of  Hermite,  it  is 
natural  to  ask  if  there  exist  weight  functions  other  than 
(1-x  )  2  for  which  a  Gauss  formula  of  Tchebycheff  type 

exists,  i.e.  which  has  equal  weights  and  which  is  exact 
for  polynomials  of  degree  <_  2n-l.  This  question  was 
answered  by  K.A.  Posse  and  N.A.  Sonin  and  the  result  is 
contained  in  the  following: 

Theorem  (Krylov  [9],  p.  184).  If  for  arbitrary  values  of 
n  =  1,2,...  there  exist  constants  A  with  x,  the  zeros 
of  the  polynomial  orthogonal  to  w(x)  on  [-1,1],  such  that 
(1.4.2)  is  exact  for  f(x)  =  1,  f(x)  =  x,  and  f(x)  =  x^  , 
then  w(x)  coincides  with  the  Tchebycheff  weight  function 
(l-x2)_5s. 

A  natural  method  of  determining  the  nodes  of  a 
Tchebycheff  quadrature  formula  is  to  set  f(x)  =  x^ 

(k=l,2 , . . . ,n)  and  the  resulting  system  of  simultaneous 
equations  will  yield  solutions  which  are  the  desired  nodes. 
Setting  w(x)  =1  in  (1.4.2),  a  direct  computation  shows 
that  the  resulting  systems  have  real  solutions  contained  in 
[-1,1]  for  n  £  7. 

1.5  EXPLICIT  FORMULAE  FOR  THE  NODES  IN  TCHEBYCHEFF  QUADRATURE 
Tchebycheff  chose  an  analytical  approach  in  determining 
the  nodes  of  the  quadrature  formula  (1.4.2).  Choosing 
f(x)  =  1,  we  are  able  to  determine  the  weight 

x  1 

A  =  —  ;  X  =  /  w(x)  dx  . 

n  -1 


. 


r 


•  i  ^  1 ::  ;  '  1 j 


Following  Tchebychef f ’ s  method,  we  choose 


f(x)  =  -  and  employ  (1.4.2)  to  yield  the  remainder 


z-x 


(1.5.1) 


R  (-M  = 
z-x 


/  dx  _  1  I  _i_ 

-1  z_x  n  k=l  z_xk 


f  SLL 


-1 


w(x)  dx  .  A  1^) 

n  n  ( z ) 


z-x 


where 


n(z)  =  (z-x1)(z-x2)  ...  (z-x  )  . 


Choosing  z  such  that  |z|>l,  we  can  write  f(x)  as 


f(x)  =  l 


x 


V 


v  =  0  z 


v  +  1 


We  can  also  write 


n ' ( z)  


l 


V 


n(z)  v=o  z 


v  +  1 


where 


n 


v 


s  =  y  xi  - 

v  k=i  k 


We  can  then  express  the  remainder  of  (1.4.2)  as 


°°  u  —  S 
In  r  Hv  n  v 


R(— )  =  l 

x  7.  -  y  L 


Z-X'  tin  __  V  + 1 

V  =  0  Z 


%  = 


/  w(x)  x  dx 

-1 


where 


’ 


(  s)  -  (c*-sM  :  -fS ) 


=  (s )  1 


3XlrlW  flBD  9W  t  X  <  I  s  I  CfjBTlcf  rfoyt  s 
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Assuming  the  quadrature  formula  to  be  valid  for  all 
polynomials  e  nn 


u  -  k  s  =  0  V  =  0,1 
v  n  v  5 

We  can  then  write  (1.5.1)  as 


jn 


# 


(1.5.2) 


w(x) 

z±x 


dx 


X  nT (z) 
n  IT  (  z ) 


00 


l 

v=n+l 


S 

V 


zv+1 


Integrating  with  respect  to  z  and  simplifying,  we  have 


°°  S  -yU  1 

(1.5.3)  n(z)  exp  [  l  — - -]=  A  exp  [y  /  w(x)  log  (z-x)  dx] 

v=n+l  v  zv  A  -1 


where  A  is  some  constant.  Since  the  expansion  of 

o  n 

00  ^  ”TV  1 

exp  [  ^  — - - -]  in  powers  of  —  is  of  the  form 

v=n+l  v  zv 

C  n  1 

1  +  —  —  +...,  and  the  expansion  of  A  exp  /  w(x)log(z-x)  dx] 

zn  -1 

is  of  the  form  A  zn  +  B  zn_1  +...,  n(z)  is  given  by  the 

polynomial  part  of  the  right  hand  side  of  (1.5.3). 


The  following  table  summarizes  the  results  of 
Tchebycheff  using  this  method  for  w(x)  =  1  and  n  <_  7 . 


-V.  0  -  V 


0  -  8  ~  -  u 


3J3  (X,.$.X)  n»lii  ftBO  9W 


1  s  :+h=v 


'  ^  '  M  /  ...„.(  X  )  v; 


(S.5.I) 


n  ‘  «  v'  :  .  x  s  oJ  3o<  ■  y  r!3.  w  ii  :  r  x x^r:  . 


.v 

r^fi 


>  nc  f.  .e.  '  >x  .2 

y-  -  5  « 

' 

■••  :  ■ 
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n 

n(x) 

nodes 

1 

X 

0 

2 

2 

X 

1 

3 

*0. 577350 

3 

X3 

X 

~  2 

0, ±0.707107 

4 

2  ,  1 

* 

±0.79^654 

4 

X 

-  3  x  +  7f5 

±0.187592 

5 

5  3,7 

r 

0, ±0.832497 

5 

X 

-  F  x  +  Y2  X 

4 

±0.374541 

±0.866247 

6 

6 

X 

_  x"  +  x2  -  1 

x  5  105 

* 

±0.422519 

±0.266635 

T 

0, ±0.883862 

7 

X7 

75.  119  3 

“  F  +  3Fo 

149 

"  FTOo  x 

4 

±0.529657 

±0.323912 

It 

was  later  shown 

that  for 

n  = 

8 

and  n  =  10 , 

the  zeros  of  the  polynomials  considered  by  Tchebycheff 
were  imaginary,  and  that  they  are  real  for  n  =  9. 


- ; i 
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1.6  IDEAS  OF  BERNSTEIN 

Tchebycheff  thus  left  open  the  problem  whether 
Tchebycheff  quadrature  is  possible  for  n  >_  10.  In 
1932,  Bernstein  proved  the  following: 

Theorem:  The  Tchebycheff  problem  cannot  be  solved  for 

n  >  9 . 

In  other  words,  the  quadrature  formula 

(1.6.1)  /  f(x)  dx  =  ^  l  f(x  ) 

-1  n  k=l  K 

does  not  hold  for  n  >  9  with  -1  <  x,  <  x0  < •  •  •  <  x  <  1. 
The  proof  results  from  the  following  lemmas  of  Bernstein: 

Lemma  1.  Let  the  formula  (1.6.1)  be  true  for  any  poly¬ 
nomial  of  degree  <  2m- 1,  m  <  n.  Then 


where  is  the  smallest  zero  of  Pm(x),  the  Legendre 

polynomial  of  degree  m. 

Lemma  2 .  If  the  Tchebycheff  formula  (1.6.1)  is  valid  for 
all  polynomials  of  degree  <_  2m- 1,  m  <  n,  then 

(1.6.2)  |  <  A<m) 

where  denote  the  Chrlstoffel  numbers  of  the  Gauss 

quadrature  formula  (1.1.2). 


, 


■ 
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■ 
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Using  the  estimates  (Natanson  [14],  pp .  522): 


JL 

m 


in  the  result  of  Lemma  2,  it  is  shown  that  for  n  >  9,  the 
inequality  (1.6.2)  fails  to  hold. 

1.7  ASYMPTOTIC  DISTRIBUTION  OF  TCHEBYCHEFF  NODES 

Although  Bernstein’s  result  established  that  not  all 
the  nodes  of  the  Tchebycheff  quadrature  formula  (1.6.1) 
will  be  real  for  n  >  9 ,  he  did  not  give  any  information  as 
to  their  location.  R.O.  Kuzmin  [10]  showed  that  if,  in 
formula  (1.5.3),  w(x)  =  1,  then  the  zeros  of  tt(x),  which 

are  the  nodes  of  the  quadrature  formula,  are  found,  with 

the  exception  of  x  =  0  for  odd  n,  in  the  neighborhood 

of  p(x)  =  1  as  n  -*  °°  where 


1 

p ( x )  =  /  log  | x-t |  dt  . 

-1 


This  result  was  rediscovered  by  Delange  [2]  who  considered 
the  location  of  the  nodes  of  the  Tchebycheff  quadrature 
formula  by  examining  the  asymptotic  expansion  of  irn(z) 


defined  by 


'  . 


-  |V  n,  ■  • 

.  I  jiri  o  j  8.1  .  )  y,  m super: I 

' 

si ' ;  c<  rus  ^  • 
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-7n  r 

Z  exp  [  -n 


Studying  the  asymptotic  distribution  of  the  zeros  of 
n  (z)  when  n  -*•  «,  it  is  shown  that  all  the  zeros,  other 
than  the  origin  for  odd  n,  are  found  arbitrarily  close  to 
a  curve  C  having  the  shape  of  an  ellipse  with  major  axis 
[-1,1]  and  minor  axis  of  length  .52552. 

1.8  TCHEBYCHEFF  QUADRATURE  ON  INFINITE  AND  SEMI-INFINITE 
INTERVAL  (Results  of  H.  Salzer) 

Although  Tchebycheff  originally  posed  his  problem 
for  the  finite  interval,  quadrature  formulae  on  the  infinite 
interval  are  equally  important.  H.E.  Salzer  [ 1 6 ]  has  con¬ 
sidered  the  Tchebychef f-Hermite  quadrature  formula 


oo 


2 


n 


(1.8.1) 


f ( x)  dx  =  A  l  f  ( x.  ) 

i=l 


and  the  Tchebychef f-Laguerre  quadrature  formula 


(1.8.2) 


/  e“x  f (x)  dx  =  A  l  f (x. ) 
o  i  =  l 


i  =  l 


For  the  Tchebychef f-Laguerre  quadrature  formula  to  be  valid 
for  polynomials  of  degree  _<  n,  it  is  necessary  and  sufficient 
to  solve  the  system  of  equations 


80*193  To  b  erl^  anlybutfS 
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n 


l  =  n( j ! )  j  =  1 , 2 , . . . ,n  . 

1  =  1  1 


The  xi  are  the  zeros  of  a  polynomial 


*n(z)  =  +  ai  zn_1  +  •••  +  an-i  z  +  an 


where  the  a^  are  defined  successively  by 


n 


n 


Jaj  +  aj-l  Xi  +  aj-2  1^1  Xi . “1  i  =  1  i  =  1 


+  ...  +  a 


n  .  ,  n 

y  x']-1  +  y  XJ 

L  l  .  L ,  l 


(j  1,2,. . . ,n)  • 

Salzer  [16]  gives  these  polynomials  explicitly  and  obtains 
their  zeros  for  n  <_  10.  It  is  seen  that  for  n=2,  x^=0 
and  X2  =  2,  and  that  for  3  £  n  <_  10,  there  exist  imaginary 
zeros . 

For  the  Tchebychef f-Hermite  quadrature  formula,  we 
are  required  to  solve  the  system  of  equations 


n 


l 

i=l 


►  ^ 

0  for  j  odd 

^  for  j  even 

2  J  /  ^ 

k  - 


J-1,2, 


,n 


As  before,  Salzer  obtains  the  nodes  by  determining  the 
polynomial  of  which  they  are  the  zeros  and  then  solves 
for  the  zeros.  For  n  =  2  we  have  x^  =  - J  ^  ,  x^  =J 


and  for  n  =  3 


nc  .  .  .  (  '  t')«  «  'i'x 


.  I-i 


. 


B  +  S 


: 
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For  3  <  n  10  there  exist  imaginary  nodes  and  hence 
Tcheby chef f-Hermite  quadrature  is  not  possible. 

Burgoyne  [1]  has  obtained  polynomials  using  the 
method  of  Salzer  for  n  <_  50  and  has  employed  Descartes’ 
rule  of  signs  to  establish  the  maximum  number  of  real  zeros 
in  the  interval  of  integration.  He  concludes  that  Tcheby- 
cheff  quadrature  of  the  type  considered  by  Salzer  is  not 
possible  for  10  <  n  _<  50. 

1.9  CONTRIBUTION  OF  GATTESCHI 

The  considerations  so  far  have  not  established  whether 
or  not  Tchebycheff  quadrature  may  be  possible  for  some  large 
values  of  n.  Bernstein’s  method,  as  outlined  earlier  in 
this  chapter,  cannot  be  applied,  since  as  pointed  out  by 
Salzer  ”...  his  method  fails  when  applied  in  exactly  the 
same  way  to  prove  the  impossibility  of  (1.8.1)  or  (1.8.2) 
(for  x^ ' s  lying  within  the  interval  of  integration)  for 
larger  n.  The  reason  is  that  it  leads  to  a  similar  in¬ 
equality,  but  now  involving  the  first  Christoffel  number 
of  the  Laguerre  or  Hermite  polynomials,  and  that  inequality 
seems  to  always  hold  for  larger  values  of  n" .  In  view  of 
the  above  remark,  it  is  interesting  to  observe  that 


■ 
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Gatteschi  [ 5 1  has  resolved  the  problem  by  a  method  very 
similar  to  that  of  Bernstein.  Denoting  by  the  constant 

of  the  Tchebycheff  formula 

b  n 

(1.9.1)  /  w(x)  f(x)  dx  =  K  l  f(x,  ) 

a  k=l 


and  C  ,  as  the  first  Christoffel  number  of  the  Gauss 
m ,  1 

formula 

b  m 

/  w(x)  f(x)  dx  =  I  c  f(xk)  3 
a  k=l  5 

Gatteschi  proves  the  following: 

Lemma.  If  (1.9.1)  has  a  degree  of  precision  2m-l,  m  <  n, 
then 

K  <  C  , 
n  m,l 

By  the  relation  (Szego  [19],  p.  360) 


C 


m ,  i 


A  ,[L'(X 
m,i  m  m 


where  A  .  are  the  zeros  of  L_(x)  . 
m,  i  m 

Similarly  in  the  Tchebychef f-Hermite  case  we  have 
(Szego  [19],  p.  360) 


J72 


2m+l 


m! 


'm,i 


[H* (h  ,)] 
m  m,i 
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where  h  .  are  the  zeros  of  H  (x) 
m,i  m 

Establishing  lower  bounds  on  l’(x  , )  and  Hf(h  ,  ), 

m  m,  1  m  m,  1  5 

Gatteschi  employes  the  preceedlng  Lemma  to  establish  in¬ 
equalities  involving  m  and  n.  The  resulting  inequalities 
fail  to  hold  for  n  >  13  in  the  Tchebychef f-Hermite  case 
and  for  n- >  9  in  the  Tchebychef f-Laguerre  case,  proving 
the  impossibility  of  quadrature  formulae  of  this  type  for 
large  n.  Combining  these  results  with  those  of  Salzer  and 
Burgoyne,  Gatteschi  concludes  that  Tchebychef f  quadrature 
is  only  possible  in  the  Tchebycheff-Hermite  case  for  n=l,2,3 
and  in  the  Tchebycheff-Laguerre  case  for  n=l,2. 

1.10  METHOD  OF  H.S.  WILF 

Wilf  [25]  has  considered  a  more  general  type  of 
quadrature  formula 

b  ,  n  b 

/  f ( x )  diKx)  =  -  l  f(x.)  ,  /  dip ( x )  =  1 
a  i  =  1  a 

Following  Wilf,  we  say  that  a  measure  dip  has  property 
T  if  Tchebycheff  quadrature  is  possible  on  a  sequence 
(nj),  (j=l,2, . . . ) ,  of  integers  tending  to  infinity,  and 
{n^ }J=1  is  called  a  T-sequence. 
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The  following  is  a  table  of  some  of  the  more  common 
measures  and  the  corresponding  T-sequences: 


di//  ( x ) 

Interval 
(a,b ) 

T-sequence 

Property  T 

dx 

ttJ  1-X2 

(-1,1) 

all  integers 

yes 

dx 

(-l.i) 

1,2, 3,  »t,  5, 6, 
7,9 

no 

e_x  dx 

(0,<») 

CM 

i — 1 

no 

2 

l —  -x 

JF  e 

dx 

(-00  ,0°) 

1,2,3 

no 

1 

l+2ax  dx 

1 — 1 

1 — 1 

1 

all  integers 

yes 

IT  J  1-X2 

1+^la^+^ax 

The  last  table  entry  is  a  result  due  to  Ullman  and  is  not  a 
commonly  used  measure  except  for  a  =  0. 

Wilf  has  also  shown  that  a  necessary  condition  for  a 
measure  to  have  property  T  can  be  determined  by  employing 
Jensen’s  Inequality,  viz: 

If  £2,...,Sn  are  non-negative,  and  if  we  define 

•  n  1 

°r  ='  I  ‘  (r=1 >2  > • • • ) 

v= 1 


then  s  <  t  implies  o  >  a. 

S  U 


n  •  •  lo  ;  z  ;  ;  i . 
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Applying  this  inequality  to 


up  =  j  xr  d iKx)  =  ^  I  x^( r=l , .  .  .  ,N) 
a  i = 1 

he  proves  the  Theorem:  Let  {n. }?  ,  be  the  T-sequence 

j  J 

of  a  measure  with  moments  (r=0 , 1 , . . . ) .  Then  for  each 
fixed  j(j=l,2,...)  the  sequence 


1 

^ p  —  ( y  2r  (r— 1j2}...j 

decreases.  If  a  >_  0 ,  then  actually 

1 

ar  =  (nj.yr)r  (r=l,2, .  .  .  ,n.  ) 

decreases . 

Applying  this  theorem  to  the  Tcheby chef f-Laguerre 
formula,  it  is  shown  that  a  increases  for  n  >  2  and 
hence  Tchebychef f-Laguerre  quadrature  is  impossible  for 
n  >  2.  By  a  similar  analysis,  Wilf  shows  that  Tchebychef f- 
Hermite  quadrature  is  impossible  for  n  _>  A ,  from  which  he 
concludes  that  the  measures  of  Hermite  and  Laguerre  do  not 
have  property  T. 

It  is  also  shown  that  if  a  measure  has  property  T, 
then  large  gaps  exist  between  the  elements  of  the  T-sequence. 
In  addition,  Wilf  conjectured  that  if  a  measure  has  property 
T,  then  it  has  zero  mass  outside  of  some  finite  interval. 
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1.11  FABER  POLYNOMIALS  AND  ULLMAN'S  CONTRIBUTION  TO  THE 
PROBLEM. 

Ullman  [23]  has  disproved  Wilf's  conjecture  by  showing 
the  existence  of  a  unit  mass  distribution  which  has  property 
T,  but  whose  mass  set  is  not  located  on  any  finite  interval. 

He  defines  a  simple  mass  distribution  of  degree  n 
to  be  a  unit  mass  distribution  with  n  equal  masses 
located  at  n  distinct  points.  The  mass  points  can  then 
be  regarded  as  the  zeros  of  the  Faber  polynomials  related 
to  the  Laurent  series  given  by 

00  mk  ai 

z  exp  (-  l  — £■)  -  z  +  a  +  —  +  .  .  . 

k=l  kzK  °  z 

and  as  the  nodes  of  a  quadrature  formula  valid  for  all 
polynomials  of  degree  <_.  n. 

Definition :  Let  f(z)  have  the  Laurent  series  expansion 

1 

ai  a2  ^ 

z  +  a  +  —  +  — p  +  .  .  .  ,  lim  sup  |ak|  <  « 
z  z  k+°° 

For  each  n  _>  1 ,  form 

[f(z)]n  =  Fn(z)  +  Rn(z) 

where  F  (z)  is  a  polynomial  of  degree  n  and  Rn(z) 
contains  only  negative  powers  of  z.  The  polynomials  F  (z) 
are  the  Faber  polynomials  associated  with  f(z).  Thus  if 

f(z)  =  z  +  fz1 2-l 


. 
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Ullman  describes  a  scheme  for  distributing  the  unit 
mass  distribution  over  the  infinite  interval  in  such  a 
manner  that  the  limiting  distribution  does  not  have  its 
mass  set  located  on  any  finite  interval  and  that  the 
Faber  polynomials  associated  with  the  distribution  have 
real  zeros.  The  nodes  of  the  quadrature  formula  are  thus 
real  and  we  have  the  existence  of  a  mass  distribution  on 
the  infinite  interval  for  which  Tchebycheff  quadrature  is 
possible . 

Ullman  [24]  has  also  shown  that  Tchebycheff  quadrature 
is  possible  on  [-1,1]  with  weight  function 


(1.11.1) 


Setting 


m,  =  /  xk  w(x)  dx 

K  -1 


Ullman  obtains  the  nodes  as  the  zeros  of  the  Faber  poly¬ 
nomials  Fn(z)  associated  with 


oo 


m 


f(z)  =  z  exp  (-  l 


k=l  kz 


. 
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For  the  weight  (1.11.1),  we  have  explicitly 


(z) 


f(z)  = 


Z+\ 


z2-l 


+  a 


z+  vz  -1 


+  a)n  +  ( 


z-\lz  -1 


n  n 


.  \  ii  li 

+  a)  -  a  . 


If  a  =  0,  then  F  (z)  becomes  the  normalized  Tchebycheff 
polynomial  of  degree  n,  and  we  obtain  Hermite’s  quadrature 
formula  whose  nodes  are  the  zeros  of  the  Tchebycheff  poly¬ 
nomial.  For  0  <_  |  a  |  ^  Ullman  shows  that  the  zeros  of 

Fn(z)  are  all  real  and  hence  Tchebycheff  quadrature  is 
possible  for  all  positive  integral  values  of  n. 


1.12  MODIFIED  TCHEBYCHEFF  QUADRATURE  FORMULAE 

We  are  thus  faced  with  the  situation  of  having 
n-point  Gauss-type  quadrature  formulae  valid  for  all 
polynomials  of  degree  <_  2n-l  and  Tchebychef f-type 
quadrature  formulae,  being  possible  for  very  limited 
n.  If  we  consider  modified  formulae  resulting  from  a 
combination  of  Gauss-type  and  Tchebychef f-type  formulae, 
two  questions  arise  naturally  and  as  formulated  by 
P.  Erdos  and  A.  Sharma  [3],  they  are: 

Problem  1:  Given  a  fixed  integer  k,  we  wish  to  determine 
n+l+k(n>k+2)  constants  ,  y^ ( i=l , 2  , . . .  ,k ) ,  x^ ( j  =  l,2 , . . .n-k) , 


and  B  so  that  the  formula 


•  '  w  <  ( 1 . 1 1 .  I  } 


/.■:  + 


'  +  ! 
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1  k  n-k 

/  f  ( x )  dx  =  l  A,  f  (y,  )  +  B  l  f  ( x  . ) 

-1  1=1  11  J  =  1  J 

is  exact  for  all  polynomials  rof  degree  <_  n+k.  We 

require  the  y1's  and  x j ’ s  to  be  in  [-1,1].  Does 

there  exist  a  number  n  such  that  for  n  >  n  the 

o  o 

formula  is  no  longer  valid? 

Problem  2:  If  for  every  n,  the  formula  of  Problem  1 

is  required  to  be  valid  for  all  polynomials  of  degree 
m  =  m(n)  <  n,  what  is  the  order  of  m(n)? 

It  is  shown  in  [3]  that  if  the  conditions  of 
Problem  1  are  satisfied,  the  answer  to  Problem  2  is 

m(n)  <  ck  [n 

where  c^.  depends  only  on  k.  Thus  Problem  1  is 
answered  in  the  affirmative  since  the  quadrature  formula 
cannot  be  valid  in  general  for  polynomials  of  degree  n+k 
if  n  is  sufficiently  large. 

A.  Meir  and  A.  Sharma  [13]  considered  the  validity 
(for  polynomials)  of  the  more  general  formula 

1  k  n-k 

(1.12.1)  /  w(x)  f(x)  dx  =  l  A  f(y  )  +  B  £  f(x.) 

-1  1=1  1=1 

where  w(x)  is  a  non-negative  weight  function.  Their 
results  are  contained  in  the  following: 


• 

• 

i’ 

• 

5 

■ 
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Theorem .  Let  k  be  a  fixed  non-negative  integer, 

2  a 

w(x)  =  ( 1-x  )  ,  a  >  -1.  Then  for  the  degree  of  exactness 
N  =  N(n,k,  )  of  a  formula  of  type  (1.12.1),  we  have 

N  <  Cn  2a+2 


where  C  =  C(k,a)  is  independent  of  n. 

Similarily  if  w(x)  =  ( 1-x ) a( 1+x ) b ,  a,b  are  real 
and  >  -1,  Jimmie  Jones  [7]  has  shown  that  if  in  (1.12.1), 
the  are  assumed  to  be  non-negative,  then  for  the  degree 

of  exactness  N  we  have 


N  £  C(k,a,b ) 


1 _ 

2 {max  ( a ,b ) }+2 


where  C(k,a,b)  is  a  constant  independent  of  n.  J.  Jones 
also  shows  that  there  exist  two  distinct  quadrature  formulae 
of  the  type 

1  3 

/  f(x)  dx  =  Af(y)  +  B  l  f(x, ) 

-1  i  =  l 


which  are  exact  for  all  polynomials  of  degree  <  5. 


.** 


CHAPTER  2 


MODIFIED  TCHEBYCHEFF  QUADRATURE  ON  (-»,«>) 

2.1  INTRODUCTION 

In  Chapter  1,  we  have  discussed  the  problem  of 
modified  Tchebycheff  quadrature  on  [-1,1]  with  weight 
function  w(x)  =  1  (Erdos  and  Sharma  [31)  and  weight 
function  w(x)  =  (1  -  x  )  ,  a  >  -1  (Meir  and  Sharma  [13]). 
However,  the  corresponding  problem  for  the  infinite 
interval  has  not  been  considered.  In  the  present  Chapter 
we  consider  the  modified  Tchebychef f-Hermite  quadrature 
formula : 


°°  2  p  n-1 

(2.1.1)  /  e-x  |  x  |  a  f(x)  dx  =  Af(y)  +  B  f(x.  ) 

i  =  l  1 

and  consider,  for  this  formula,  problems  of  the  type  posed 
in  Section  1.12  for  the  finite  interval.  Our  results  can 
be  formulated  by  means  of: 

Theorem  1.  If  a  formula  of  type  (2.1.1)  is  valid  for  all 
polynomials  of  degree  <_  2m+l,  m  <  n,  with  real  A,  B,  y 
and  x^(i=l,2, . . . ,n-l) ,  and  a  >  -h  ,  then 

m  <  C  log  n 


where  C  =  C(a)  is  independent  of  n. 
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The  case  when  a  =  0  is  of  some  interest  in  itself, 
being  connected  with  Tchebycheff  quadrature  using  the 
common  Hermite  weight  function  on  the  infinite  interval. 

We  then  state: 

Corollary  1.  Assuming  the  formula 

°°  2  n-1 

/  e  f (x)  dx  =  Af (y )  +  B  £  f(x. ) 

i  =  l  1 

to  be  valid  for  all  polynomials  of  degree  <_  2m+l,  m  <  n, 
with  real  A,  B,  y,  and  x^ ( i=l , 2 , . . . ,n-l ) ,  then 

m  <  C  log  n 

where  C  is  a  constant  independent  of  n. 


2.2  PRELIMINARIES 

In  order  to  prove  Theorem  1,  we  shall  introduce  the 
polynomials  (x)  which  are  a  natural  generalization 

of  Hermite  polynomials  (Szego  [19]  P.  377).  H^a^(x)  are 

orthogonal  polynomials  corresponding  to  the  weight  function 
e~x  |  x  |  a,  a  >-h  in  [-°°,°°].  Their  relation  to  the  Laguerre 
polynomials  is  given  by 


(2.2.1) 


( a ) 


m  ~2m 


H2m  (x)  =  (-1)U‘  2  m!  Lm 


(a-h)  ,2 


(x  )  ,  a 


T r  ( ot )  /  \  ,  ,  Nm  02m+l  ,  T  (  a  +h)  ,  2  >.  , 

H2m+1  x  =  2  m  Xm  ^x  >  >  a  >  ~%  • 


\ 
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Denoting  the  zeros  of  H^a\x)  by 


^  <  4*1  <*••<  4*1 

1  s  m  2  ,m  m  ,m 


(  a)  ,  \ 

and  the  zeros  of  L  (x)  by 

m  J 


<  „<«)  <...<  u(») 

l,m  2,m  m,m 


we  have  from  (2.2.1)  and  Szego’s  estimate ([ 19 ]  p.  118) 

( a ) 

on  the  largest  zero  of  the  Laguerre  polynomial  1/  (x) 


j 


(2.2.2) 


>  ( 2m+a-3/2 ) 


'2m  ,2m- 


( a) 


m,m 


^2mh,2m+D2  =  ^  >  <2ra+“-^  ' 


m,m 


V 


Combining  the  estimates  of  (2.2.2)  we  obtain 


(2.2.3) 


<  -  (m+a-3/2)^ 

l,jn 


We  shall  also  require  Stirling’s  formula  for  the 
Gamma  function: 


r(m+l)  <  2  J  urn  (^-)m  . 


(2.2.4) 


r 


■ 


1 


n:  t  n 


( . ;) 


1 


(£\£-»+fn)  -  >  3 
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2.3  SOME  LEMMAS 

We  shall  require  the  following  lemmas  in  the  proof 
of  Theorem  1. 

Lemma  1.  If  (2.1.1)  is  valid  for  polynomials  of  degree 
_<  N,  then  B  >  0 . 

Proof.  Examining  (2.1.1)  we  see  that  N  >  2  and  choosing 


f(x)  = 

X 

1 

'C 

a 

\s 

we  have 

°o  2  0 

n-1 

0  < 

r  -x  |  1 2a  /  n  2 

J  e  | x |  (x-y)  dx  =  B 

I 

— 

00 

i  =  l 

Since 

(x,-y)2 

>  0 ,  we  must  have  B  >  0. 

Lemma  2 .  If  formula  (2.2.1)  is  valid  with  A,  B,  y  and 
x^ (i=l,2, . . . ,n-l)  for  polynomials  of  degree  <_  N,  then 
it  is  valid  for  polynomials  of  degree  <_  N  with  A,  B, 
-y  and  -x^ (1=1 ,2 , . . . ,n-l)  . 

Proof.  Choosing  f(x)  =  xv  ,  0  <_  v  <_  N,  we  have 

(2.3*1)  /  e”x  | x | 2a  xv  dx  =  Ayv  +  B  \  xV  . 

i=l 

For  v  an  even  integer,  the  Lemma  is  obvious.  If  v 

00  2  2a  v 

is  an  odd  integer  we  have  /  e~x  |x|  a  xv  dx  =  0  so  we 

oo 

can  replace  y  by  -y  and  x . (i=l , 2 , . . . ,n-l )  by 
-x^(i=l,2, . . . ,n-l) .  Lemma  2  is  thus  proven. 


•  :  .1  ' 


( y  ) 

- 
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(  a ) 

Lemma  3 .  Let  (x),  a  >  -  Jg,  denote  the  orthogonal 

polynomials  of  §2.2  and  let  ,  be  the  zeros  of 

- — -  -  i,m  i  =  l - 

H(ot)(x)  . 
m 

If  y  >  0  and  R(x)  is  the  polynomial  of  degree 


2m+l  given  by 


(2.3.2) 


R(x)  = 


[H^a)(x)]2(x-y)2 

(x-£U) ) 

^x  n,m; 


m  >  2 


then 


(2.3.3) 


00  2  p 

/  e_X  | x |  a  R(x)  dx  <  0  . 


Proof.  R(x)  may  be  written  as 

(2.3.4)  R(x)  =  [H^)(x)]2(x+4“h2y)  +  (c[“ty)2  ^ 

k  ^l,m' 

2 

Since  [H^a^(x)]^  and  e"*x  |x|2a  are  both  even,  it 
m  11 

follows  that 


(2.3.5) 


x  |  x  |  2ot  x[H(a)(x)]2  dx  =  0 

ii  m 

—  oo 


I  e' 


By  the  orthogonality  of  the  H^a^(x)  , 


(2.3.6) 


/  e 

.  oo 


2  ,  [H^a) (x)  ]2 

-x  i__i2a  m  _ 

(x-chb 

l,m 


x 


dx  =  0  . 


:  i  -  /  _ 


. 


I  y - —  -  (x)H  (S;£.S) 
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Multiplying  (2.3*4)  by  e  x  |x|^a  and  integrating, 
we  obtain 

OO  2  QQ  2 

(2.3.7)  /  e-x  |x|2“  R(x)  dx  =  (t,\a)-2y)  /  e_x  |  x  1 2o[H^“ )  ( x)  ] 

-00  J-5m  _oo  171 

(  n  ) 

For  m  >  2,  £1  '  <  0,  whence  the  result  (2.2.4)  follows. 

—  ’  l,m 

Lemma  4 .  If  formula  (2.1.1)  is  exact  for  all  polynomials 
of  degree  <_  2m+l,  m  <  n,  and  if  the  abscissae  are  given 
by  x]_  <  x2  <***<  xn  i>  then  we  have 

( a ) 

x,  <  Cn 
1  1  ,m 

Proof.  By  Lemma  2,  we  may  assume  without  loss  of  generality 
that  y  >_  0 .  Then  for  the  polynomial  R(x)  defined  by 
(2.3*2)  we  have  from  (2.3*3) 

00  2  p  n-1 

(2.3*8)  0  >  /  e_x  |xpa  R( x )  dx  =  B  £  R(x.)  . 

i  =  l 

Therefore,  since  B  >  0  by  Lemma  1,  there  exists  a 

v  such  that  R(x  )  <  0.  From  (2.3*2),  R(x)  <  0  only  if 

x  <  *  Hence  x  <  for  some  v  and  x,  <  • 

l,m  v  l,m  1  l,m 

2.4  PROOF  OF  THEOREM  1 

Proceeding  with  the  proof  of  Theorem  1,  we  can  choose 


dx  . 


y  >  0  without  loss  of  generality  by  Lemma  2.  We  will 


. 


•  ‘do  evr 
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consider  the  proof  In  two  cases  determined  by  the  location 
of  y . 


Case  1.  Suppose  |  x-^y  |  >_  lxn-1-y|  . 


x — v 

Choosing  f(x)  =  ( - *-)  and  using  (2.1.1)  we  have 


xi_y 


n-l  x . -y  , 

(2.4.1)  nB  >  B  l  (-i— r  =  /  e 

i=l  xl"y 


-x2.  | 2a  ,  x-y  N  2 


dx 

x1-y 


00  -X 

I 


1  x 1 2a(x2+y2 ) 

(x-^y  )2 


00  -t ,  a .  .  2  N  ,, 

r  e  t  (t+y  )  dt 

J  .  .  p 


(x]_-y) 


dx 


4  : 

;  > 

"I'V 


r(a+3/2)  +  y ^  r(a+^) 

(xx-y)2 

C^d+y2) 


(xfy) 


From  (2.4.1),  we  obtain  |  a  lower  bound  on 

I  ' 


b  : 


(2.4.2) 


B  > 


C, (l+y*) 

_ ± _  • 

n(x^-y ) 2 


We  now  choose  f(x)  =  ( :;X‘“^-) 2  ( ~) 2m  ^  and  again  use 
formula  (2.1.1)  to  obtain 


xi_y  xi 


. 

-  K 

4 

: 

. 


9  \  -  (—=)  a  <  Bn 

*  -  ~if  -± 


.  '  $ — 


( y~ x) 


■ 

(SlC+I)j? 
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(2. It. 3)  /  e_x‘|x|2a  (-*ZZ)2 (JL)2m~6  dx  =  B  nf  (^ii.)2(^i) 


xi-y 


i=l 


xl"y 


x_. 

x. 


2m- 6 


>  B  . 


Proceeding  to  evaluate  the  left  hand  side  of  (2.4.3) 
we  obtain  the  following  upper  bound  for  B  in  terms  of  m! 


7  e-x2|x|2°  (^rZ)2(JL)2m-6 

l  1  1  xi-y  xi 


0°  —  x^  |  |  2a  /  2m-4,  2  2m-6>. 

dx  =  r  e _ L*] _ (a  .  1  ±y_x _ 1  dx 

ax  j  “  “2  2nP5  ax 

(x1-y)  xx 


(2.4.4) 


r ( m+a- 3/2 ) +y  r(m+a-5/2) 

<  - T -  2  -2"m“6 

(x  -y)  x 


c_2 ( 1+y  ;  r (m+a-3/2 ) 

72.  2m- 6 

(xx-y)  x± 


Employing  (2.2.3)  and  Lemma  4,  we  have 

(2.4.5)  x^  <  _  (m+a-3/2 )^ 

Using  Stirling’s  formula  for  the  Gamma  function  and 

(2.4.5) ,  we  can  express  (2.4.4)  as 


. 
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)  dSM  '-  'o  \ 


at  — 


(  \  -v+m)  .  ■  y.+  (sv:-o+m )i  _ 

;  x  ■ 


'*  .  fi  J 


1  y-  >•  1 
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•  <• 
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7  -x2  i  i2a,  x-yN2,  Xs2m-6  ,  ^3 

J  e  |  x  |  (- — £■)  (— )  dx  <  — 


C„(l+y2)(m+a-5/2)m+a’2 


xi-y  X]_ 


m+a-5/2,  ,  n/oNm-3/  \2 

e  ^  (m+a-3/2)  (x^y) 


(2.4.6) 


Clt(l+y2)(m+a-5/2)a+1 

m3  72 

e  (x1-y) 


Combining  (2.4.2),  (2.4.3)  and  (2.4.6)  we  have  the 
following  inequality  involving  B: 


(2.4.7) 


( 1+y2 ) (m+a-5/2) 0+1  C  (1+y2) 

-  >  B  >  - o 


m,  ,  2 
e  (x1-y) 


n(x1-y ) 


Simplifying  (2.4.7)  we  have 


(2.4.8) 


Cc (m+a-5/2 ) 
2 _ 

m 


a+1 


1 

>  — 
n 


Taking  logarithms  of  (2.4.8)  yields 


Cr  +  (a+1)  •  log  (m+a-5/2)  -  m  >  -  log  n 


which  can  be  simplified  to 


(2.4.9) 


m  <  log  n 


where  =  C^(a)  . 


■ 
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Case  2.  Suppose 


x 


n- 


1 


-y 


> 


X 


1 


-y 


Since  y  >  0,  we  must  have 


lx  .  1 

> 

1  x,  | 

> 

|  ha)  1 

H 

1 

C 

1  1 1 

1  1  ,m  1 

Proceeding  as  in  Case  1,  we  choose  f(x) 
arrive  at 


x-y 


xn-l-y 


and 


(2.4.10) 


B 


cg(i+y2) 

n(xn.i-y) 


2 


If  we  then  choose 


f(x) 


x-y  n  2 ,  x 


xn-l-y 


)  ( 


X 


n-1 


2m-6 


we  obtain,  as  in  Case  1,  the  following  inequality: 


(2.4.11) 


Cg ( 1+y 2 ) (m+a-5/2 ) a+1 

m  /  n  2 

e  xn-l_y 


>  B 


Combining  (2.4.10)  and  (2.4.11)  and  simplifying,  we 


have 


(2.4.12) 


C10 (m+a-5/2) 0+1 


m  n 

e 


•  !&-!  <  !;*'  <  lr-,rv-1 

■ 

c  •'  f. 

(OI,t».£) 
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Taking  logarithms  of  (2.4.12),  we  obtain 

(2.4.13)  m  <  log  n 

where  C  =  C-^(a)  . 

Choosing  C  to  be  the  larger  of  in  (2.4.9) 

and  C  in  (2.4.13)  completes  the  proof  of  Theorem  1. 


sie.'iw 
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CHAPTER  3 


QUADRATURE  FORMULAE  FOR  THE  SEMI-INFINITE  AND  INFINITE  INTERVAL 

3.1  INTRODUCTION 

In  this  chapter,  we  will  consider  the  infinite  Tchebycheff 
quadrature  formula 

°°  2  0  n 

(3.1.1)  /  e”x  | x |  a  f (x)  dx  =  B  l  f  (x, ) 

i  =  1 

and  the  semi-infinite  Tchebycheff  quadrature  formula 

00  n 

(3.1.2)  /  e”x  xa  f(x)  dx  =f  B  l  f(x.  )  ♦ 

o  i  =  1  1 

Wilf  [25]  and  Gatteschi  [5]  have  considered  the  case 
a  =  0,  and  their  results  are  summarized  in  Chapter  1.  However, 
they  have  not  considered  the  degree  of  exactness  N  of  these 
formulae.  Our  estimates  for  N  are  contained  in  Theorems  1 
and  2,  which  give  a  qualitative  generalization  of  the  results 
of  Gatteschi  and  Wilf  for  the  infinite  and  semi-infinite 
intervals . 

3.2  INFINITE  INTERVAL 

Theorem  1.  If  a  formula  of  type  (3.1.1)  is  valid  for  poly¬ 
nomials  of  degree  <_  2m- 1,  m  <  n,  with  real  B  and  real 
x^(i=l,2, . . . ,n) ,  and  a  >  -^ 


then 
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m  <  K  log  n 


for  some  constant  K  =  K(a)  . 

We  will  require  the  following  lemma  for  the  proof  of 
Theorem  1  and  will  use  the  notation  of  §2.2. 

Lemma  1 .  If  formula  (3.1.1)  Is  valid  for  polynomials  of 
degree  <_  2m- 1,  and  If  the  ( 1  =  1 , 2  , .  .  .  ,n)  satisfy 


x 


1 


<  X 


2 


<  •  •  •  <  x 


n 


then 


Proof.  Choosing  f(x)  = -  -  ,  a  polynomial  of  degree 


[H(a)(x)]2 

m 


(  a ) 

2m- 1,  we  can  employ  the  orthogonal  property  of  the  Ir  (x) 
and  formula  (3.1.1)  to  obtain 


(3.2.1)  0  =  /  e"x 


00 


2 


X 


Since  [H^a^(x.)]^>  0  and  H^a^(x.)  4  0,  for  all  i 
ml—  mi'  s 


( m )  (a) 

(x^-S^  ^)  <  0  for  some  xj_  •  Hence,  (x^-^  <  0 


yields  the  desired  inequality 


which 


r  "  [  X  >  m 


1:   ‘K 

' 
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Proof  of  Theorem  1. 

Returning  to  the  proof  of  Theorem  1,  we  will  now 

r  (  q  \ 

establish  an  upper  bound  on  the  constant  —  — ~ -  in  terms 
of  m. 

Since  we  have  assumed  (3.1.1)  to  be  exact  for  poly- 

x  2  m  2 

nomials  of  degree  <  2m-l,  we  may  choose  f(x)  =  (— ;  ~  . 

xi 

Then  from  (3.1.1),  we  have 


(3.2.2) 


/  e“x  | x | 


—  oo 


XI 2c(jL)2m-2  dx  .  T(a±h\ 

X1  n  i=l 


n  x.  0  0 

l  (^)2m-2 


rU+ki . 


n 


Evaluating  the  left  hand  side  of  (3.2.2)  we  have 


(3.2.3)  /  e 


—k  i  i 2o ,  x  s 2ra-2 
X  . 

X, 

-oo  1 


oo  -t,m+a-3/2 

/As  c\\\  —  ^  ,  ret  , , 

(— )  d  x  =  / - - - - -  dt 


o 


2m- 2 


r  (m+ot-3§)  . 

2m- 2 

x. 


Employing  Stirling's  formula  in  (3.2.3)  we  obtain 


(3.2.4) 


2  a 


2m- 2 


dx 


.  o/o \m+a-l 
C-L(m+a-3/2 ) 

<  m+a-3/2  2m-2 

e  x1 


Combining  the  estimates  of  (3.2.2)  and  (3.2.4)  yields 


V. 
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C, (m+a-3/2 )m+a  ^  r(a+^) 

(3.2.5)  — — r — 5-75- — -5--  y  ■  >  -  # 

em+a-3/2  x2m-2  n 

Employing  the  estimate  (2.2.3)  and  Lemma  1,  we  have 


(3.2.6) 


-(m+a-3/2),s. 


Substituting  the  above  inequality  in  (3.2.5)  yields 


(3.2.7) 


r ( a+h) 


(m+a-3/2 ) 


m+a*1! 


n 


!m+a-3/2(m+a-3/2)ra-1 


C, (m+a-3/2 ) a ' 

<  -  * 

m 

e 

Taking  logarithms  of  (3.2.7)  yields 

C^i  -  log  n  <  Cp.  +  (( a:)  )  log  (m+a-3/2)  -  m 

which  can  be  written  as 

(3.2.8)  m  <  K  log  n 

where  K  =  K(a),  completing  the  proof  of  Theorem  1. 


3.3  SEMI-INFINITE  INTERVAL 

•  Theorem  2 .  If  formula  (3.1.2)  is  valid  for  polynomials  of 
degree  <_  2m- 1,  m  <  n,  with  real  B  and  with  real  positive 
x^(i=l,2, . . . ,n) , a  >  -1  then 


•  --7-  < 
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m  <  K  log  n 

for  some  constant  K  =  K(a). 

We  shall  require  the  following  lemma  in  the  proof  of 
Theorem  2  and  will  use  the  notation  of  §2.2. 

Lemma  1 .  If  formula  (3.1.2)  is  valid  for  all  polynomials 
of  degree  <_  2m- 1,  m  <  n,  and  if  the  nodes  x.  (i=l,2 , .  .  .  ,n) 
satisfy 


X-,  <  x0  <  •  •  •  <  x 

12  n 


and  a  >  -1,  then 


(3.3.D 


X  >  y 
n 


(  a ) 
m,m 


[La)(x)]2 

Proof.  Choosing  f(x)  =  — - — 7 — \ ,  a  polynomial  of 

(x-Lah 

m  ,m 

degree  2m- 1,  we  can  employ  the  orthogonal  property  of  the 

L^a\x)  and  formula  (3.1.2)  to  obtain 
m 


ia)(xi)]2 


00 

r ( a+1 ) 

n 

[L 

r  _ —  X  a 
J  e  x 

f(x) 

dx 

- 

l 

0 

n 

i=l 

( 

(X^]2 

1  0 

and 

[L(a)(x. 
m  l 

)]  + 

0 

(  a) 


l  m,m 


( a ) 


( a ) 

(x.  -  y  )  >0  for  some  x.  . 
l  m,m  i 


(  a ) 

Hence  (x  -y  )  >  0  which  yields  the  desired  inequality 
n  m3m  J  H  J 


(a) 

x  >  y 
n  m  ,m 


•2oot5 


Xt>  (x)l  °XX  9  \ 


0  <  [(.xV’M] 
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Proof  of  Theorem  2. 

We  now  proceed  with  the  proof  of  Theorem  2  using  the 
method  of  Theorem  1. 

Since  (3.1.2)  is  assumed  to  be  exact  for  polynomials 

of  degree  £2m-l,  we  can  choose  f(x)  =  (— j  ~  and 

n 

obtain  from  (3.1.2) 

(3.3.2)  /  e-xx“(i)2m-1  dx  =  l 

o  xn  n  i=l  xn  n 

We  also  have 

r ( 2m+a )  . 

2m- 1 

xn 

Using  Stirling’s  formula  in  (3.3*3),  the  resulting 
inequality  is 

,  0  ,  -i  \  2m+a-^ 

1(2m+a-l) 

2m+ a - 1  2m- 1 

x 

n 

Combining  Lemma  1  and  (2.2.2)  we  have 

x^  >  2m+a-l. 

Substituting  this  inequality  in  (3*3.4)  yields 


(3.3.4) 


r  -x  a  /  x  >,  2m-l  , 

/  e  x  (— )  dx 

o  n 


(3.3.3) 


r  -x  a  /  x  >,  2m-l  , 

J  e  x  (-— )  dx 

A 

o  n 


-  -f.'  -•  ;  v- 
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00  0  C0(2m+ot-l) 

/  0  0  r- \  r  -x  a,  Xx2m-1  ,  2 

(3.3.5)  J  e  x  (— )  dx  < 

A 

o  n 


2m+a-^ 


2m+a-l / 0  ,  , v  2m- 1 

e  (2m+a-l) 


C3(2m+a-l) 


a+^ 


2m 


Combining  (3*3.2)  and  (3*3.5)  we  have 


(3.3.6) 


C1|(2m+a-l) 

2m 


a+h 


>  i  • 

n 


Taking  logarithms  of  (3.3.6)  gives  the  inequality 


(3.3.7) 


m  <  K  log  n 


where  K  =  K(a)  which  completes  the  proof  of  Theorem  2. 


* 
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CHAPTER  4 


PARTIAL  RESULT  AND  OPEN  PROBLEM 

4.1  INTRODUCTION 

In  this  chapter,  we  consider  the  modified  Tchebycheff 
quadrature  formula  on  the  semi-infinite  interval  (0,»): 

oo  n-1 

(4.1.1)  /  e"x  f ( x )  dx  =  A  f(y)  +  B  £  f(x.) 

o  i=l  1 

and  propose  the  same  problem  as  for  the  corresponding 
formula  (2.1.1)  for  the  interval  (-«>,«>)  in  Chapter  2. 

Our  results  in  this  connection  are  partial  and  incomplete. 
Although  we  believe  that  the  analogue  of  Theorem  1, 

Chapter  2  holds  true  for  formula  (4.1.1),  we  have  not 
been  able  to  prove  this  so  far.  Hence  the  analogue  of 
Theorem  1,  Chapter  2  is  formulated  here  as  a  conjecture 
in  Theorem  2. 

One  of  the  reasons  for  the  failure  of  an  extension 

of  our  method  of  proof  of  Theorem  1  of  Chapter  2  is  that 

we  can  no  longer  determine  the  location  of  y  in  relation 

to  the  larger  nodes  of  the  above  quadrature  formula.  In 

xi-y 

other  words  we  cannot  establish  that  ( - )  <  K  where 

K  is  a  constant  independent  of  n,  for  some  x^  sufficiently 
large,  and  hence  we  can  not  establish  inequalities  of  the 
type  (2.4.4). 
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We  shall  now  prove  the  following: 

Theorem  1.  If  (4.1.1)  is  exact  for  all  polynomials  of 
degree  <_  2m+l,  m  <  n,  with  A  and  B  real,  y  and 

xi-y  p 

x4 (i=l , 2 , . . . ,n-l)  contained  in  (0,«),  and  (- - — )  <  C 

1  1  1  X  -i  "V 

n-1 

where  C>1  is  a  constant  independent  of  n,  then 

m  <  K  log  n 

for  some  constant  K  independent  of  n. 

4.2  SOME  LEMMAS 

We  shall  require  the  following  lemmas  in  the  proof 
of  Theorem  1. 

Lemma  1 .  If  (4.1.1)  is  valid  for  polynomials  of  degree 
<_  N ,  then  B  >  0 . 

Proof.  Examining  (4.1.1),  we  see  that  N  >_  2  and  choosing 

2 

f(x)  =  (x-y)  gives  the  result. 

Lemma  2 .  The  degree  of  exactness  N  of  formula  (4.1.1) 
is  a  non-decreasing  function  of  n . 

Proof.  Suppose  formula  (4.1.1)  is  valid  for  polynomials 
of  degree  *  <  N.  Then  defining 


A(n+1) 

=  A(n) 

-  B(n)  , 

B(n+1)  _  B(n) 

(n+1) 
x : 

=  x<n) 

( l<i <n-l) 

l 

1 

(n+1) 

(n) 

(n+1) 

X 

n 

=  y 

,  y 

y 

... 
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we  see  that  a  formula  of  type  (4.1.1)  with  n  replaced 
by  (n+1)  is  valid  for  all  polynomials  of  degree  <_  N 
and  hence  the  degree  of  exactness  is  a  non-decreasing 
function  of  n. 

Let  the  zeros  of  the  Laguerre  polynomial  L  (x) 
be  denoted  by 


<  y . 


l,m  2,m 


<  •  •  •  <  y 


m,m 


and  let  L  (x)  be  normalized  such  that 
m 


9 

/  e”x  L  (x)  dx  =  1  * 
J  m 

o 


Lemma  3 .  If  (4.1.1)  is  exact  for  polynomials  of  degree 


v 


m  ,m 


<  2m+l,  and  y  <  — then 

_  2  - 


X  n  >  y 

n-1  m,m 


Proof.  Consider  the  function  R(x)  = 


From  (4.1.1)  we  have 


L2(x) (x-y ) 2 
m 

(x-y  ) 
m  ,m 


(4.2.1) 


/  e  x  R(x)  dx  =  B  l  R(x. )  . 

i=l  1 


n-1 


o 


Since  we  can  write  R(x)  as 


R(x)  =  L  (x)[x-y+(y  -y ) ( 1+ 

m  J  m,m  J  x-y 


y  y  -y 

m,m  )  y(Mm,m 


m,m 


x-y 


)] 


m  ,m 


M  >  93  Tg  '  -  1  r  'l  ■ 

.  ...  . .  ■  •• 

h  cton*.  ■ 06 

n t  rn  • 

j  rl  rioui  .•  '  t  .Hon  9J 
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we  have 


/  e“  R(x)  dx 
o 


7  e-xL2(x)[x-y+(ll  -y )  ( ('^CL)  ] 

j  m  J  pm,m  J  x-y  J  x-y 

o  5  m,m  m,m 


(4.2.2) 


uu  _  oo 

/  e_xL2 (x )  xdx  -  2y  /  e”xL2(x)  dx 


o 


m 


+  y _  /  e”xL2(x)  dx 


m  ,m 


o 


>  -  2y  +  y 


m,m 


>  0  . 


Hence  we  have  from  (4.2.1) 


n-1 


(4.2.3) 


B  l  R( x .  )  >  0  . 
1  =  1  1 


Since  B  >  0  by  Lemma  1,  (4.2.3)  will  be  true  only  if 


(4.2.4) 


n-1  L2(x. ) (x.-y)2 

l  -S—i - i -  >  0 

i=l  (x.-y  ) 

i  m,m 


and  (4.2.4)  will  hold  only  if  x  >  yv.  _  for  some  v  and 
hence 


v  ■  m,m 


X  -»  >  y  4 

n-1  m,m 


■  •  .  ■  '  . 


'  '  \ 


(S.S.t  ) 


u  <+  v*,^  -  .< 


•  ■  :  : 
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Lemma  4 .  If  (4.1.1)  Is  exact  for  all  polynomials  of  degree 
<  2m+l,  5  m  <  n,  then 


xn-l  > 


m  ,m 


Proof.  If  y  <  — we  can  use  Lemma  2  to  obtain  the  desired 

2 

inequality . 

ym  m 

We  will  now  consider  the  proof  for  y  >.  — *—  .  Let  us 

y  2 

assume  that  xn_]_  ±  •  We  have  (Szego  [191  pp.  127) 

y  <  4m+3  and  (Szego  [19]  pp.  126) 

111  y  ill 

tt2  (m-5/4 ) 2  , 
ym-l,m  4  in+1/2 


Since 


4m+3  tt  2  (m-5/4)^ 


m+1/2 


m  >  5 


we  have 


m  .m 


<  y 


m-l,m  ' 


Choosing  f(x)  = 


Lm(x) (x”y) 


,  a  polynomial  of  degree 


2m- 1,  we  have  from  (4.1.1) 


(4.2.5) 

oo 

0  =  / 


e"XL^(x) (x-y ) 


dx  = 


O  (x-y _ )  (  X-y 


m,rrr 


m-1  ,m 


) 


n-1 

B  l 


L  ( x.  ) ( x . -y ) 
m  i  l  J 


i=l  ( x . -y  ) (x-y  ) 

l  m ,m  I  m-1 ,m 


I! aol  dsfixe 

i.>  n;  "  ' 


beiiB9b  sriJ  riicddo  od  S  BmmbJ  esu  nso  9W  "-L  .  >  ^  II  .'ioon<I  . 

• 

-  >•  ^ 
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The  right  hand  side  contains  only  negative  terms  which  is 
impossible.  Hence  we  have  a  contradiction  and 


4.3  PROOF  OF  THEOREM  1 

By  Lemma  2,  it  is  sufficient  to  prove  the  theorem 
when  m  >  5  and  we  choose  n  in  (4.1.1)  to  be  a  sufficiently 


large  integer.  Choosing  f(x) 


have  from  (4.1.1) 


c1d-y+y2) 


and  hence 


(4.3.2) 


C2(l-y+y2) 


Choosing  f(x)  =  ( - —  )2( - 

&  x  .,-y  x 

n-1  J  n-1 

in  (4.1.1)  we  have 


-)m  ^  and  substituting 


' 


orftf  sritf  9voTq  otf  Snetoll'iUB  al  .£  smmeJ  yd 
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We  now  obtain  an  upper  bound  for  B  by  finding  an  upper 
bound  for  the  left  hand  side  of  (4.3.3)  as  follows: 


/  e~x(— '-)c(^~) 


x-y 


m-3 


x  -y 
n-1  J 


x 


n-1 


d  x  =  / 

o 


°°  -X/  m-1  m-2  .  2  m- 3 n 

e  (x  -yx  +y  x 


(x  •, -y )  (x  -,)m  ^ 
n-1  J  n-1 


dx 


(4.3.4) 


C3  r (m) ( 1-y+y  ) 

,  >.2/  \  m—  3 

(xn-l-y)  (Vl> 


Combining  Lemma  4  and  the  bound  on  y^  m  given  by  (2.2.2) 
we  have 


(4.3.5) 


1 


2m- 1 
>  -  • 

4 


Using  Stirling’s  formula  for  the  Gamma  function  and  (4.3.5) 
in  (4.3.4)  we  have 


(4.3.6)  /  e“A( 
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”X' — — )m-3  dx  < 


x  -i  —  y  x 

n-1  J  n-1 


C4(m-l)m-^(l-y+y2) 

2  m-1 , (m-h) sm-3 


(xn-l-y)  e  ( 


) 


Combining  (4.3.2),  (4.3.3),  and  (4.3.6)  we  have 


Cj-  (m-1)  5/2  ( 1-y+y2  )2m  3  Cp(l-y+y2) 

9  -  >  B  >  -6 - *■ 


(  \  2  m 

xn-l-y  e 


n(xn-l-y) 


which  can  be  simplified  to 


(4.3.7) 


Cg  m5/2  2m-3  x 

m  >  n 
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Taking  logarithms  of  (4.3-7)  yields 

C  ^  +  5/2  log  m  +  (m-3)  log  2  -  m  >  -  log  n 

which  can  be  written  as 

m  <  Cg  log  n  • 

This  completes  the  proof  of  Theorem  1. 

4.4  OPEN  PROBLEM 

On  the  basis  of  Theorem  1  and  the  corresponding  Theorem  1 
of  Chapter  2,  we  formulate  the  following  unproven  theorem: 
Theorem  2.  (conjecture) 

If  (4.1.1)is  exact  for  all  polynomials  of  degree 
_<  2m- 1,  m  <  n,  with  A  and  B  real,  and  y  and 
x. (i=l,2, . . . ,n-l)  in  (0,«),  then 

m  <  K  log  n 

for  some  constant  K  independent  of  n. 

This  conjecture  is  also  offered  when  the  weight  function 
in  (4.1.1)  is  e~x  xa  9  a  >  -1. 
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